Analytic Continuation of the Doubly-periodic Barnes Zeta Function 

Guglielmo Fucci*' and Klaus Kirsten^^ 

^Department of Mathematics, East Carolina University, Greenville, NC 27858 USA 

^Department of Mathematics, Baylor University, Waco, TX 76798 USA 

(Dated: April 16, 2013) 

CO , The aim of this work is to study the analytic continuation of the doubly -periodic Barnes zeta func- 

^^ ■ tion. By using a suitable complex integral representation as a starting point we find the meromorphic 

(N ■ 

extension of the doubly periodic Barnes zeta function to the entire complex plane in terms of a real 



< 









integral containing the Hurwitz zeta function and the first Jacobi theta function. These allow us to 
explicitly give expressions for the derivative at all non-positive integer points. 

I. INTRODUCTION 

The Barnes zeta function, introduced for the first time in ||li |2|], represents a higher dimensional gener- 
alization of the Hurwitz zeta function 



^H{s,a) = Yj{n + ay'' for5l-s>l. (1.1) 

^ . Namely, let 5' e C, /i e R+, and r e R.^. For %s > d the Barnes zeta function is defined through the series 

^, 

^' : ^s(s,fi\r) - y Cu + n • rr , (1-2) 

o ' w 



and it can be analytically continued in a unique way to a meromorphic function in the entire complex plane 
possessing only simple poles at s = 1,2, ...,d. In this work we shall be interested in the meromorphic 
j^ , extension of a function closely related to the two dimensional Barnes zeta function. Let b,c eM.+, Q.m,n = 

ibm + en with m,n eZ, and a e C\n„, „. For ''Ks > 2 we consider the following zeta function 

({s, a, b,c) = 2, C*^ + i^'" + cny" , (1.3) 

(m,«)EZ2 

which is the kind of zeta function resulting from Dirac operators on the two-torus as considered in gener- 
alized Thirring models [8]. The zeta function defined above is doubly-periodic with respect to the variable 
a. In fact for /: € Z one has the relations ^(s, a + kc, b, c) = ^(s, a, b, c) and ^{s, a + ikb, b, c) = ^{s, a, b, c), 
and due to this double-periodicity we assume, without loss of generality, that < %a < c and < 3a < b. 
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Since the zeta function (11.31) is analytic in the semi -plane ^^^ > 2 and the ratio of the two periods c and \h 
is not a real number, one can conclude that (ll.3l l belongs to the class of elliptic functions [11]. 

The main idea of the present work is to represent the doubly-periodic Barnes zeta function (ll.3l l for 
2^5^ > 2 in terms of a contour integral in the complex plane. The desired analytic continuation of (11.31 ) to 
the region 51 5 < 2 is then achieved by a suitable deformation of the integration contour. This process yields 
an expression for (ll.3l l valid in the entire complex plane in terms of an integral over the interval [0, 1] of 
the Hurwitz zeta function and the logarithmic derivative of the first Jacobi theta function. The analytically 
continued expression for ^{s,a,b,c) will allow us to very easily compute its values at all integer points, 
s eZ. In addition, we will also provide an explicit expression for the derivative of ^(s, a, b, c) with respect 
to s at all non-positive integer points. 

We would like to point out that one of the main advantages of our study of the analytic continuation of 
(11.31) is that its double-periodicity property remains manifest in all the formulae. This is an aesthetically 
pleasing feature that can also be desirable if one wishes to implement these expressions in a computer pro- 
gram for numerical evaluations. The study of a zeta function closely related to the doubly-periodic Barnes 
zeta function considered here has appeared in ^ where a method for obtaining its analytic continuation 
was used which differs from the one we employ in this work. 

The outline of the paper is as follows. In the next section we construct a contour integral representation 
for ^(s, a, b, c) valid for %s > 2. From this representation we perform, in Section In] the analytic contin- 
uation in the entire complex plane. Section JV] is devoted to the explicit computation of the derivative of 
({s, a, b, c) with respect to the first variable at all negative integer points and at i' = 0. 

II. CONTOUR INTEGRAL REPRESENTATION OF ^(s, a, b, c) 

As we have already mentioned in the Introduction, the main idea of our approach is to represent the 
doubly-periodic Barnes zeta function (11.31 ) in terms of a contour integral. Since by assumption c ?^ we 
introduce, for convenience, the function 

fmin, s) = , (2.1) 

(a„, -I- ny 



with n eZ and a,,, = - + i-m € C\{0}, where m e Z. Obviously, in terms of the newly introduced function 
fmin, s), the zeta function (11.31 ) reads 

(is,a,b,c)^c~'^^fyn{n,s) . (2.2) 



By utilizing Cauchy's residue theorem we rewrite the sum over the index n in (12.21 ) as a contour integral. 
More precisely one has 

y, 7 = Yjmin, s)^ :r-r \ Az fmiz, s) ncot{nz) , (2.3) 

where F is a contour that encloses counterclockwise all the (simple) poles of the function 7rcot(;rz). Let us 
point out that the representation (I2.31 l is well defined in the region ''Ks > 2. 

Before specifying the integration contour Y in detail, we would like to observe that the function fmiz, s), 
obtained from (12.11 ) by replacing n with z, possesses branch cuts extending from the points z = -am- 
The exact position of the cut will depend explicitly on the summation index m. First, note that from the 
assumptions stated below (11.31) one obtains the inequalities b/c > 0, %{a/c) > 0, and 3{a/c) > 0. This 
allows us to conclude that ^^m > for all m, da^ > for m > and 3 am < for m < 0. The last remark 
shows that the cut lies in the lower half complex plane when m > and in the upper half complex plane 
when m < - 1 . The contour F has to be chosen in such a way to enclose only the poles of cot(7rz) but not 
the branch points of the function fmiz, s). More precisely, the contour is the union F = F+ U F_ where F+ 
satisfies the property < 3F+ < -So-m for m < -I while F_ satisfies 3a„, < 3F_ < for m > (which 
simply means that the contour is closer to the real axis than the cut). The contour F is depicted in Figure [T] 
with -a>o and -a<_i denoting representatives of the set of branch points with index m > and m < -1, 
respectively. 

With the contour of integration completely determined we can express the left hand side of (12.31) as a 
sum of two contributions 

y 7 = T-^ I dz fniz, s) ncot(nz) + t-^ I dz fm(z, s) ncot(nz) . (2.4) 

^ (am + ny 2m Jr+ 2;ri Jy_ 

The next step of our approach is to rewrite the integrand in (12.41 ) in a way that will allow a separate 
treatment of the integral over F+ and over F_. By utilizing the complex exponential representation of the 
function cot(7rz) and after some straightforward algebraic manipulations we obtain 

cot(7rz) - -i f 1 + ^^- 1 , for 3z > , (2.5) 



e-2mz _ 1 

and 



cot(7rz) = i 1 1 + -T- I , for 3z < . (2.6) 



glmz _ I 

Since 3F+ > and 3F_ < 0, we use the representation (12.51 ) for the integral over F+ and the representation 
(12.61 ) for the integral over F_. Proceeding in this fashion leads to the results 

^ J dz fmiz, s) jTcotinz) - -^ J dz fmiz, s)- j dz ^^^, (2.7) 
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Figure 1 : Integration contour used for the representation (12.3 



and 



2m Jy_ 



dz fmiz, s) n coXijiz) 



- I dz /„,(z, s) + I 



fmiz, S) 



(2.8) 



The next step is the defomiation of the integration contours which is the subject of the following Lemma. 

Lemma 1. Let T = Y+^T^be the integration contour described above, and z = u + a^ with m eZ. Then, 
in the semi-plane ^s > 2, one obtains 



^ J^dz Uz, s) ncotinz) = - [du ^^^ _ ^ + f 



Au 



g2m(u-a>o) _ J ' 



(2.9) 



where y is a contour enclosing in the clockwise direction the negative real axis including the point u = 0. 

Proof. The proof of this result is based on a suitable deformation of the integration contour F. Before 
deforming the contour, we focus on the first term on the left hand side of (12.71) and (12.81 ). From the definition 
of the function fniz, s) in (12.11 ) we can write 



I dz fmiz, s) - I dx ix + ani + i3F+) ■* . 



(2.10) 



Since, by assumption, the analysis is restricted to the region %s > 2 and since 2la,„ > one can conclude 
that for all m 



i: 



dz fniz, s)^0 



(2.11) 



By using a similar argument one can prove that for all m the following relation holds 

dzf,n{z,s)= dx (x + a„, - i3r_r^ = . (2.12) 

Let us consider next the second integral on the right hand side of (12.7b . which can be rewritten as 



g2jTBr 



^— Y [(X + a>())2 + (3r+)2] ^ ^ , (2.14) 



r fmiz, s) ^ r ^^^^^ i3r+)-^^ — i— - — , (2.13) 

and is convergent for ''Ks > 2. If the branch cut extends from the points -a>o, which simply means that it 
lies in the lower half plane, the contour r+ can be shifted away to infinity in the upper half plane and does 
not contribute. In fact 

(x + «>() + i3r+)~^' : —= 

as 3r+ — > oo. If, on the other hand, the branch cut extends from the points -a<-i in the upper half plane, 
we shift the contour r+ around the branch cut as shown in Figure |2] 

The second integral over r_ in (12.81 ) is convergent for 51 i' > 2 and can be expressed as 

r dz 4^ - r dx (x + a^- idr.y——-^—- . (2.15) 

By using arguments similar to the ones outlined above one can prove that if the cut extends from the branch 
points -Q'<_i in the upper half plane, then we can shift the contour r_ away to infinity, namely 3r_ -^ -oo, 
and the integral (12.151 ) vanishes. If, instead, the branch cut extends from the points -a>o, in the lower half 
plane, we shift the contour r_ as shown in Figure |2] 

We can therefore conclude that for m < -I only the integral over the deformed r+ gives a non- vanishing 
contribution and by making the substitution u = z + a<-i we obtain 

where y is a contour enclosing in the clockwise direction the negative real axis including the point m = as 
shown in Figure |3] By using a similar argument, when m > 0, we obtain 

2^ I dz Uz, s) n coti^z) = j^ du ^-^—^^ . (2. 17) 

By substituting the results (12.161 ) and (12.171 ) in the relation (12.41) we find the claim ( 12.91 ). n 

Lemma 1 immediately allows to rewrite the doubly-periodic Barnes zeta function. First note, that from 
the expressions (12.11 ) and (12.21 ) we can write 

C,{s,a,h,c) = —-\ \ dzfm{s,z)cQMnz) , (2.18) 
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Figure 2: 
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Figure 3: 



which by using Lemma[T]gives 



^{s, a, b, c) - —c 



-oo _ 

m=-i ^y 



-2m(H-f-i|m) 



oo _ 

1 m=o ^y 



,2m(M-f-i|m) _ J 






'y 



g2m(»-f ) _ J 



+ c 



oo ^ 



du u " 



1 



1 



(2.19) 



_j ^y ,^g2m(«-f)+2;r|m _ ^ g-2;:i(«-f )+2;:|m _ J J ' 

where the integral in the second line represents the contribution due to m = 0. By introducing a new variable 
^ e R defined as ^ = e~'^^^'^ and by combining the two exponentials in the last line of (12.191) we find the 
expression 



^(s,a,b,c) = c 



I 

1 ^y 



,2m{u—al c) _ 1 



Au u 



-2i sin 2n 



('"h^Di 



q 



2m 



2cos(27r[M- f])(?2'« + ^4m 



, (2.20) 



which, once again, is valid for ''Ks > 2. The second integral appearing in (I2.20l i. although rather cumber- 
some, can be expressed in terms of a simple special function. In fact, let t € C with \t\ < 1 and z e C; the 
first Jacobi theta function 7?i(z; t) has the following representation as an infinite product Ig] 



7?i(z; t) - 2t^/4 sinz ]~[ (l - 2t^" cos(2z) + /") (l - r^") . 



(2.21) 



n=[ 



By taking the derivative of the logarithm of (I2.21l i and by using q introduced below ( 12.191 ) for r, we obtain 
the formula ||5|] 

\ni'}<(7:a^ = ^^^'^'' - cot 7 + 4s.m(27^ ^^ ■ — - 

dz 



d &[(z;q) ^ 

— In i?i(z; q) = ^\ ^ ^ cotz + 4 sin(2z) 2_j : 



(2.22) 



&i{z;q) " ' ~'f^ 1 - 2^2n cos(2z) + ^'*" ' 

Since |^| < 1, for any finite z the sum in ( 12. 22b is absolutely convergent |19|]. This imphes, in particular, that 
the convergence of the series in (12.221 ) is uniform in z in the region |3z| < nb/c. The last remark justifies 
the interchange of the sum and the integral in (I2.20l i to obtain, by using the expression (12.221 ) with the 
substitution z — > niu - a/c), the representation 



^(s,a,b,c) = c 



2m{u~alc) 



I 

--— c \ du u -— ln??i TT 
2n Jy du \ 



1 2 Jy 



du u "^cotlTr 



a 

u 

c 



a 

u 

c 



;? 



(2.23) 



where the contour y is chosen so that |9y| < nb/c in order to allow the use of (|2.22| i. By combining the first 
two integrals in (12.231) and by using the relation 

1 



i 1 

+ - cot TTX = , 

e2mx _ 1 2 2 



(2.24) 



we find 



l{s,a,b,c) = -^r-c ' \ du u -— In i?i tt m ;q\ — — \ du u . (2.25) 

2n Jy du \ y c\ I 2 L 



Since %s > 2, the second integral in (12.251 ) gives a vanishing contribution and we are left with the following 
compact representation 



^{s,a,b,c) = ---c ■' I du u "-—Ini^iin u ;q] 

2k Jy du \ y c\ J 



(2.26) 



We would like to point out that the last integral representation preserves the double-periodicity property of 
the original sum in (11.31 ). In fact, let m, n eZ.By using (12.261) we write 



({s,a + nc + imb,b,c) = c ^ \ du u ' — Inj^ilTr 

27T Jy du \ 



a 

u 

c 



nn - min—\q 



(2.27) 



Now, let Zm, n = Qn + nT)n denote the vertices of the fundamental parallelogram in the z-plane. Then the 



first Jacobi theta function is quasi-periodic on the lattice 

&^{z + (m + nT)n-T) = (-l)'"+"e-'(2"^+™'"^7?i(z;T) . 



(2.28) 



By utihzing (IHSl) in (IZTTl) we obtain 

^(s, a + nc + imb, b, c) = —nc'" I du u~^ + ^(s, a, b, c) = ^(s, a, b, c) , 

Jy 



(2.29) 



where, as before, the last equality follows from the fact that since %s > 2 the integral over y vanishes. 

The representation (12.261 ) allows to easily compute ((s, a, b, c) for s = n, n e N, n > 2. For these values 
of s, the contour encloses a singularity at m = of order n and the residue theorem immediately shows 

d 



^{n,a,b,c) c " \ du u " — lni?i(7r 

2;t Jc^ du \ 



a 

u 

c 



;^ 



1 d" 



cl {n- ly.du" 
n\" 1 d" 



ln'&i{u;q) 



h\-&i{u;q) 



.cl {n- ly.du" 
We next use ( 12.261 ) to construct the representation of ({s, a, b, c) valid in the whole complex plane. 



(2.30) 



III. ANALYTIC CONTINUATION 



The integration contour y in ( 12.261 ) consists of a union of three distinct paths, namely -y = y+ U Ce U y_, 
where Cg is the circular portion of radius e centered at the origin, y+ is the straight path positioned at a 
distance 5 above the negative real axis, and y_ represents the straight path positioned at a distance 5 below 



9 



the negative real axis. Furthermore, for later use, we denote by -e the projection on the negative real axis 
of the intersection of 7+ (or y_) with the circular portion Ce of the integration path. This remark allows us 
to rewrite (I2.261 l as 



({s,a,b,c) = --T-c ] du u *-— ln&i{n u ',q]- ■;^c * I 

2n Jc^ du \ y c\ / 2n Jy^uy_ 



du U ■ -—In 7?i I TT 

du 



a 

u 

c 



\q 



(3.1) 



The last representation is a suitable starting point from which we can proceed with the analytic continuation 
to the region %s <1. 

The first term in (13. IK namely the integral along C^, is left unchanged; it will be dealt with later. The 
second term in (13. Il l can be expressed as a sum ly^ + ly , where ly^ represents the integral over the path 
7+ and ly the one over the path 7_. We will present details for the integral ly^, ly follows accordingly. 
We parameterize the path y+ as m = i(5 + x, with x e (-00, -s\, and rewrite the integral over y+ as follows 



ly. = 



= c ^ I du u ■' — In??! In u ;q] 

In Jy^ du \ y ci / 

2n J_o. 



dx (x + 16) ' — ln??i \n 
dx 



X + is 

c 



dx(i(5-xr' — ln??i TT 



is - X 

c 



(3.2) 



27r Jg ^ ' dx 

The integral in the last line of (13.21 ) is well defined at the lower limit of integration. To establish the 
convergence at the upper limit of integration it seems we would need to analyze the asymptotic behavior of 
the integrand as a; — > 00. However, this analysis can be avoided by using the following quasi-periodicity 
property 



■&i{u + n;q) = -j9'i(m; q) , 



(3.3) 



which is obtained by setting n = and m = I in the more general formula (12.281 1. 
In order to exploit the above property we represent ly^ as 

°° rn+l+E J 



Iy.= — c~^ > I dx {is - xy^ — Xn-^iin is - X •,q\ 

In ±riJn+s dx \ y ci J 



n=0 

and perform the change of variables x = y + nto obtain 

— c-^y dyiiS-y-nr-ln&An 

In ^,Jg dy \ 



I 



7+ 



is — y — n 



(3.4) 



(3.5) 



Next we send 5 ^ 0, in which case e — > £. By using the quasi-periodicity (13.31 ) in the expression (13.51 ). we 
then have 



00 /-.!+£ 1 



n=0 ' 



dy 



a 

-y-- 
C 



;<? 



(3.6) 
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We note, at this point, that the series Z^gCy + n) ' converges unifomily in y in the interval [s, 1 + e] for 
51^ > 1. This allows us to rewrite (13.61) as 



I 



7+ 






l+£ 



dy ^H{s,y)—lnd^i\n 



a 

-y-- 
C 



;^ > 



(3.7) 



where the Hurwitz zeta function (h(^, x) in equation (11.11) has been used. As is well known, the Hurwitz 
zeta function can be analytically continued in a unique way to a meromorphic function in the entire complex 
plane possessing only a simple pole with residue 1 at the point s = I. 

For the integral J^ over the lower part of the contour we proceed similarly to obtain 



J. 



7- 



In Je 



1+^ d 

dy ^H{s;y)—ln&i{7r 



a 

-y-- 
C 



By adding the contribution from ly^ and ly we obtain the result 



-* y+Vy- ~ ■'■ y+ +-'-}'_ — 



sin(7r5) 



l+e 



dy l^H'yS\y)—\rv&\\n 



a 

3;+- 

c 



(3.8) 



(3.9) 



Noticing that the integral along C^ in ( 13.11) is defined for i^ e C, the above expression performs the analytic 
continuation of l^{s, a, b, c) in terms of the Hurwitz zeta function ^h{s, y) to the full complex plane and it is 
valid for £ > small enough. By choosing %s < 1, the limit e ^ can be performed and the integral along 
Cs in ( 13.11 ) can be seen to vanish. 

We can now summarize the results obtained in this Section as follows: 

Theorem 1. Let ({s, a, b, c) denote the doubly-periodic Barnes zeta function defined in di.il ). 
If%s < 1, then 



sin(7r5') _.. r^ d / 

i^{s,a,b,c) = c' Jy ^//(5;y)— ln??i TT 

^ Jo Ay \ 



a 

y+- 
c 



\q\ ■ 



If s = n with n e N+ — { 1 ), then 



(in, a,b,c) -\ — 



1 d" 



■lni^i(M;^) 



(3.10) 



(3.11) 



yd {n- l)!dM" 

Due to the prefactor sin(;ri')/7r in the integral representation (I3.10l l. the substitution s = -n with « € No 
leads to a vanishing contribution of ^{s, a, b, c) when the first argument is a negative integer including zero. 
This remark proves the following 

Corollary 1. Let n e No, then ^(-n, a, b, c) = 0. 

In order to find the representation of (,{s,a,b,c) valid for 2ls > 1, note that for j ^ we have the 
behavior 



^H{s,y) = -^+U{s) + 0(y), 

y' 



(3.12) 



11 



where the first term results from the n = contribution in (ll.ll ). This is the term responsible for the 
restriction %s < 1 found in the previous results. The analytic continuation of (13.101 ) to the full plane is 
found by observing that 



Ms,y)^ —+^His,y+l) 



(3.13) 



By substituting (13.131) in the representation (13.101 1 we find that the second term in (13.131 ) gives a function 
holomorphic for s e C The resulting contribution to (13.101 1 coming from the first term in (13.131 1 gives, 
instead, the integral 

1 



= I dy y '—Ind-iin 



a 

y + - 

c 



;^ . 



(3.14) 



which is valid for 51 i^ < 1. At this point, partial integration can be applied repeatedly to obtain a representa- 
tion that extends expression ( 13.141 ) further to the right of 51 5 = 1. For example, after one partial integration 
(I3T41 ) reads 

1 



/n = 



1 



l-s 



d 

dy 



(^3' + ^;'?) - j dyy '^^ — In&iln y + ^ ■,qj 



(3.15) 



which is now valid ior ''Ks < 2. After n partial integrations a representation valid for 51 5 <« + 1 is found 
and (13.111 ) can be verified from there. The above representation also shows 



n d 
^(l,a,b,c) = --—ln&i{u;q) 
c du 






(3.16) 



where the first term comes from /q, and the second term from the pole of ({s;y + 1) at 5^ = 1. The second 
term can be simplified noting 



In ; ^ = —17T , 



'^iH;-?) 

which can be proved starting from the property (13. 3t . 



(3.17) 



IV. THE DERIVATIVE OF THE DOUBLY-PERIODIC BARNES ZETA FUNCTION AT NEGATIVE 

INTEGERS 



In addition to the results already presented, from the integral representation (I3.101 l one can compute 
the derivative of the doubly -periodic Barnes zeta function with respect to the first variable. Differentiating 



12 



a 

y+- 

c 



(13.101 ) leads to the result 

t{s,a,b,c) - cos{ns)c~^ I (iv ^//(5;v)-— In i^i I tt 

Jo dj V 

sin(7ri') _.. r' d / 

- c Inc (i)' ^//(s;)')-— In j^i TT 

n Jo dj V 

d 

d^ 



a 

y + - 

c 



sin(7ri') _, 
+ c 



7T 



dy CH(^;y)— In &i\n 



a 

y + - 

c 






(4.1) 



where, here and in the rest of this work, the prime indicates differentiation with respect to the variable s. By 
setting s = in ( 14.11 ) and by noting that for n € No the following relation holds ia] 

B„+i(x) 



^H{-n,x) ^ — 



n + I 



(4.2) 



with B„(x) denoting the Bernoulli polynomials defined in terms of the Bernoulli numbers S^ as 

B„ 
we find 



ax) = J](^)B,x"-^ 



(4.3) 



C{-n,a,b,c) = {-\y 



"■4;("r);^X'*/--^-.(" 



a 

y+- 
c 



(4.4) 



Integrating by parts (14.41 ) yields the result 



(_l)«+ic" 

C\-n,a,b,c) = ln??i(7r 

n+ I 



a 

- + 1 
c 



k=Q 



1 \^k + -. Dn+\ in 



(-1)"+V' 



n + 



r" Z (" r)^" ~ ^ ^ 1)B^ jl' dj /-^ In ,^1 (;r 



3^ + - 
c 



;^ ■ 



(4.5) 



At this point it is convenient to distinguish between two cases: n - and n > I. The reason for this 
distinction lies in the relation ||6|] 



zf-rw 



for « > 1 



(4.6) 



k=o \ ■■ ' I Bo for n = 



By exploiting the above relation and (13.171 ). one obtains 



^'(0,a,&,c) = -ln?9'i (tt-;^) + — + I dy In&iin y + - ■,q] 



(4.7) 



and, when « > 1 , 



(-l)"c" 



« + 1 



(-l)"c" f^/« + l 



a-n, a, b, c) = in'-^^Bn+i + '-^ J] ( " ^ ' )(« - /c + 1)B^ | dy y"''' In &i | 



« + 



k=Q 



\n-k+\)Bk f dyy"-Hn&i(ny+- ■,q\ 



■q] .(4.8) 



The integrals that appear in (I4.7l i and (I4.8l i can be computed from the results of the following lemma. 
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Lemma 2. For < 3z < 7t3t with q = e'"'^, one has a Fourier-type representation 



1 V — 1 CI cos(2?2/) 

\n'9\{z;q) - -\nq + Iht/Ct) + ln(2sinz) - 2 ) 

o ^ I - q 



,2n 



(4.9) 



Proof. The logarithmic derivative of the first Jacobi theta function can be represented in terms of an infinite 
series for < 3z < ttSt as 13,160 

-2" — q^" sin(2z) 

Anti-differentiation of ( I4.10I ) yields 

q^" cos{2nz) 



— lni^i(z;fl) - cotz-i-4 y r- sin(2«z) - cotz + 4 V - — 

dz ^ \ - q^" '^ y - 

n=\ ■* n=\ 



Iq^"" cos(2z) + (f"" 



(4.10) 



In 79-1(2;^) = /(^)+lnsinz-2 



y Ji__ 

n=\ ■■ 



- /(^)+ln sin z+Y_^ In (l - 2^^" cos(2z) + ^'*") , (4. 11) 



n=l 



where /(^) is an arbitrary function. In order to determine the unknown f{q) we use the infinite product 
representation ||5l|6|] 



T?i(z; q) = 2G{q)q^ sinz f] (l - 2^^" cos(2z) + /") , 



(4.12) 



H=l 



where 



G(^)-]~[(l-^2„J (4^3) 

From equation (I4.12l i and by recalling that the Dedekind eta function is defined for Or > and q = e'"'^ as 

00 
rj{T) = e'^Y\{^-q^'') , (4.14) 

n=l 

one can easily find 

lni?i(z;^) - ln2 + ln77(T) + - ln<7 + Insinz + 2ln(l - 2q^" cos(2z) + ^^") . (4.15) 

By comparing (14.111 ) with (14.151 ) we finally find that 

(4.16) 
Substitution of the explicit expression (14.161 ) for the function f{q) in (14.111 ) yields the claim ( 14.91 ). n 



/(^) = ln2 + ln?7(T) + -ln^. 



The use of ( 14.91 ) in the expression (14.71 ) allows us to obtain the following result for the derivative of the 
doubly-periodic Barnes zeta function at 5 = 0, 

^'{0,a,b,c) = -ln&i(7:-;q\ + — + -lnq + lnr]ii-\+ I Inllsin nly +-]] . (4.17) 
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When « > 1, we use once again the result (|4.91 l in (I4.8l l. By perfomiing this substitution and by recalling 
4.61) one obtains 



^'{-n,a,b,c) = 



(-l)"c" v^/« + l 



tI 



(-l)"c" „ 
+ 1^ — fi«+i, 



(« - A: + \)Bk 



OO '^/ 1 



« + 1 

where we have introduced for typographical convenience the functions 

-1 



„_, (^) = r d). y"-' In (2 sin ["()■+ ^)]) 



dy y" '' cos (inl 



a 

y + - 



(4.18) 



(4.19) 



(4.20) 



Note that the integral appearing in (14. 17| i reduces, according to the definition ( 14.191 ). to Iq (a/c). 

The integrals (14.191 ) and (|4.20| i can actually be explicitly computed in terms of polylogarithmic functions 
and trigonometric functions, respectively. The following lemma provides an expression for the function 

In-k {ale). 

Lemma 3. For « e No and for A e C\Z, 



In{A) 



+ 



l{n + l){n + 2) n+l 



> I K — I )■ , 



n + 



n+\\, _i\-^ k\ / i 



|; (^ - ')a A/,,, (e^-^) . (4.21) 



l(k-l 



Proof. Performing an integration by parts leads to the expression 

J„ (A) - ^^ + ^^ + In sin(7rA) — \ Ay /+' cot {n \y + A]) . (4.22) 

n+l n+ \ n+ \ n+ \ Jq 

The integral in (14.221 ) containing the cotangent can be explicitly computed in terms of polylogarithmic 

functions. First, by exploiting the change of variables j + A — > j, we obtain 

jl dj f*' cot {n [y + A]) - j; (" ^ ^)(-A)"+i-^ £''' dy / cot {ny) , (4.23) 

and by rewriting the cotangent in terms of complex exponentials the integral on the right-hand side of (I4.23t 
takes the form 

.2mv 



I 



■A+1 i^A+l i^A+\ 

dy / cot (ny) = -i ( dy y^ + 2i \ 
Ja Ja 



dyy' 



e2my _ I ■ 

This can be rewritten using the polylogarithmic function defined for \z\ < 1 and s £ C by the sum 



Li.a) = Z ^ ' 



n=l 



(4.24) 



(4.25) 
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and by analytic continuation in the entire complex z-plane |l7|]. Of particular importance to our analysis is 



the following property satisfied by the polylogarithmic function 



H 



^Au„p") = u„,(.-). 



and since 



Uo(e^'''y) = 



„lmy 



glniy _ I ' 



we have that 



Jlmy 



_J_ALii(.2->). ^. 
Ini Ay ^ ' e^my _ | 



The result (14.281 ) employed in (14.241 ) yields 






(4.26) 



(4.27) 



(4.28) 



(4.29) 



Integrating by parts k times and using, at each step, the relation (I4.261 l leads to 



-r 



■^*'d„>Au,(.-).-i2^ 



1 ^ k\ (-i)"y 



n.,k—n 



«=0 

k 



(k-n)\ {Inx) 



i^Lwi (.-'>■) 



A+l 



In light of the previous result and after computing the elementary integral on the right-hand side of (I4.24t 
we can conclude that 



-A+l 



r d. / cot (..) = -^ E n 'V""' --y.^Tr^yx~. "V'-'u.^ [e^-^] 



(4.31) 



By substituting expression (14.311 ) in (14.231 ) we obtain, for (14.221 ). 



In 2 in 1 i7r(-A)"+2 "^ /„ + i\(-i/+i v^ A + 1\ / 

/„(A) - - + - + -lnsin(7rA) + ^-^— y i^-^ V A"' 

« + 1 w+l « + 1 n+l ^\ A: /Ar+l -^ \ / / 

^^§("r)<-4^(^)'(|rrH-(-)- 



(4.32) 



The above result can be simplified further, in fact 

m{-AT*^ y /« + 1\(-1)^ y Ik + 1\ _/ _ i7r(-A)"+2 y (n + \\{-\f^' 



k=0 



1=1 



and by using the relation, valid for a € 



m, 



k=o 



k k+l 



k+\ 
l.-l -1 



, (4.33) 



y lm\^^ _ {a + 1)'"+' 
2-i[kjk+l ~ m+ 1 



(4.34) 
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one obtains 



'^^tirwtr^'^ 



m 



(n + l){n + 1) 



(4.35) 



*:=0 ^ ' l=\ 

The I - contribution from the last term in (I4.32l i can be simplified by following an argument similar to 

ll 



the one leading from (I4.33l l to (I4.35I ). By noticing that 



Lii(z)--ln(l-z), 



(4.36) 



one obtains 



{-Af 



n «+i 



n+l 



mA In 2 



1 37ri 

In sin(7rA) - . (4.37) 



n+l /-J\ k P 4i\j) ^ ' n+l n+l n+l ' ' 2{n + I) 

k—\J j—i 

By substituting (14.351 ) and (14.371) in the expression (14.321 ) the claim (14.211 ) immediately follows. 



D 



Let us now focus on the analysis of the function j7„-/t,/(a/c) in (I4.20I ). which can be rewritten, after a 
suitable change of variables, as 



^"-^•'(?) = ?d)^'^°^( 



cos I Inl- 



Inl 



I 



dx X cos x + 



sin 1 2nl- \ \ dxx'^ sin x . 



(4.38) 



(2;r/)*+i 
The use of known reduction formulae for the trigonometric integrals on the right-hand side of (14.381 ) yields 



/ 



dxx 



/ 



cosx 



sm;K: 



\ 



k\ 



y ^^(-1)' 



X*="2' + 



COSX 



[¥] 



k\ 



y — - — (-1)' 



k-ii-\ 



smx J ^—^ (k -21 - 1)! 



, (4.39) 



I sin X J \ - cos X j -^ {k - 21) ! 
where [x] denotes the integer part of x. The application of the expression ( 14.391 ) to (14.381 ) provides the result 



' \c/ (27r0*^+i \ cl ^{k-2j-l)r \ 2 ) 



{2nlf 
I 



-ljcT-\k\ 



{2niy 



1+1 



sin 1 2-nl- 



V— ^— (-l)^(2;r0*^-2;+/ 
^(^-27)!' I 



l+(-l)* 



-ii%\ 



(4.40) 



We can conclude that the equation (14.181 ). together with the explicit expressions (14.211 ) and (14.401 ). gives 
a formula for computing the derivative of the doubly-periodic Barnes zeta function at all negative integer 
points. Moreover, an expression for ^'{0,a,b,c) is obtained from ( 14.171 ) by using ( 14.211 ) with n = 0. The 
results obtained in this section allow for a very efficient way of computing ^'{-n,a,b,c) as the explicit 
formulae can be easily implemented in an algebraic computer program. For completeness, we display the 
expression for ^'{-n, a, b, c) when n = 0,n = -I, and n = -2. 

For n = 0, one has 

^'(0,a,b,c) = -ln&i(7T-;q\+-lnq + lnr]ii-\+mi---\, (4.41) 
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for « = -1, we obtain 

"^ -2' 1 / a^ 



and for n = -2, one gets 

2 2 °^ 21 ^ 

a-2,a,b,c) = -^^U,(e'-"^)-'-Y,j^pCOs(2nl^) . (4.43) 

The above expressions are valid in the range of parameters stated below equation (11.31 ) and have to be 
periodically continued, as is clear from the fact that the integral representation (|4.4| i is periodic. 



V. CONCLUDING REMARKS 

In this work we have performed a detailed analysis of the meromoiphic extension of the doubly -periodic 
Barnes zeta function in the entire complex plane. The contour integral representation (I2.3l l allowed us, 
after a suitable deformation of the integration contour, to obtain the meromorphic extension of the doubly- 
periodic Barnes zeta function ^{s, a, b, c). This analytically continued expression revealed to be particularly 
appropriate for the explicit computation of the values ({n, a, b, c) and of the derivative ^'{-n, a, b, c) for 
n e No. 

The process of analytic continuation delineated in this work is quite general and its applicability is not 
limited exclusively to the study of the doubly -periodic Barnes zeta function. In fact, the method developed 
here can be applied to more general elliptic functions. Although the representation of elliptic functions in 
terms of integrals over R+ has been constructed, for instance in ||3|], our method would provide a contour 
integral representation for the class of elliptic functions. Since we have seen that such representation has 
some advantages, such as the double-periodicity of the results and the almost straightforward computation 
of the values and derivative at specific points, it could, perhaps, provide either new results or simplify 
already known ones regarding elliptic functions. This seems to be an interesting idea which deserves further 
investigation. 

Aside from their intrinsic mathematical interest, our results can find applications in problems related 

o 

to physical systems. The Thirring model is used to describe simple interacting field theories HlOll . The 
one-loop partition function for generalized Thirring models is proportional to the derivative at s = of the 
doubly-periodic Barnes zeta function ||8|]. The result (14.411 ) obtained here can then be directly applied to the 
analysis of these interacting field models. 
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